function [ s ] = evaluateQuadraturePoint( p, q )
%EVALUATEQUADRATUREPOINT Summary of this function goes here
%   Detailed explanation goes here

global poissonRatio shearModulus;

shiftAlpha = 1e-6;

r = p - q;

rn = max(norm(r), 1e-6);

rg = r / rn;

M = [rg(1)*rg(1), rg(1)*rg(2); rg(2)*rg(1), rg(2)*rg(2)];

I = diag(ones(2, 1));

s = -1 / (8 * (1 - poissonRatio) * shearModulus * pi) * ((3 - 4 * poissonRatio) * log(rn + shiftAlpha) * I - M);

end

